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We study different architectures for a photonic crystal in the microwave regime based on super- 
conducting transmission Unes interrupted by Josephson junctions, both in one and two dimensions. 
A study of the scattering properties of a single junction in the line shows that the junction be- 
haves as a perfect mirror when the photon frequency matches the Josephson plasma frequency. We 
generalize our calculations to periodic arrangements of junctions, demonstrating that they can be 
used for tunable band engineering, forming what we call a quantum circuit crystal. Two applica- 
tions are discussed in detail. In a two-dimensional structure we demonstrate the phenomenon of 
negative refraction. We finish by studying the creation of stationary entanglement between two 
superconducting qubits interacting through a disordered media. 
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I. INTRODUCTION 

Circuit QED-^ is quantum optics on a superconduct- 
ing chip: a solid state analogue of cavity QED in which 
superconducting resonators and qubits act as optical 
cavities and artificial atoms. After successfully repro- 
ducing many key experiments from the visible regime 
— qub it-photon strong coupling and Rabi oscillations^, 
Wigner function reconstruction^, cavity- mediated qubit- 
qubit coupling^, quantum algorithms^ or Bell inequali- 
ties measurement^ — , and improving the quality factors 
of qubits and cavities, c-QED establishes as an alterna- 
tive to standard quantum optical setups. 

The next challenge in the field is the development of 
quantum microwave photonics in the Gigahertz regime. 
The scope is the generation, control and measurement 
of propagating photons, contemplating all its possibili- 
ties as carriers of quantum information and mediators 
of long distance correlations. The natural framework is 
that of active and passive quantum metamaterials^ with 
open transmission lines to support propagation of pho- 
tons and embedded circuits to control them*"^^. Qubits 
can be a possible ingredient in these metamaterials. A 
two-level system may act as a saturable mirror for res- 
onant photons^^^^, as it has been demonstrated in a 
breakthrough experiment with flux qubits^^, continued 
by further demonstrations of single photon transistors^^, 
and electromagnetically induced transparency^^. These 
groundbreaking developments, together with theoretical 
studies of band engineering using qubits*"^ and Joseph- 
son junction arrays^ ^, and recent developments in the 
field of photodetection^^^^, provide solid foundations for 



this rapidly growing field. It is important to contrast 
these developments with alternative setups in the high- 
energy microwave regime (Terahertz) which differ 
both in the architecture and the scope. 

In this work, we advocate an alternative architecture 
for both passive and active quantum metamaterials based 
on transmission lines with embedded Josephson junctions 
(JJ). Adopting a bottom- up approach, we first study the 
scattering of travelling photons through a single junction, 
the simplest and most fundamental element in supercon- 
ducting technologies. It is shown that, in the few photon 
limit, the linearized junction acts as a perfect mirror for 
resonant photon. Starting from the single JJ scattering 
matrix, we show how to engineer metamaterials using 
periodic arrangements of junctions both in one and two- 
dimensional transmission line networks. Compared to 
previous approaches, this combines the travelling nature 
and flexible geometry of photons in transmission lines^^, 
and instead of qubits^'^'^^ it relies on the simple and ro- 
bust dynamics of a linearized junction^^. Previous pro- 
posals lacked one of these two ingredients. 

The simplicity of this setup opens the door to multi- 
ple short-term applications. In this manuscript we dis- 
cuss mainly two. The first one is the observation of a 
negative index of refraction in a two-dimensional circuit 
crystal. This would be achieved by injecting an appropri- 
ate microwave in a square network of transmission lines, 
where only half of it is populated with embedded junc- 
tions. Second and most important, we study the interac- 
tion between qubits in a disordered quantum metamate- 
rial, showing that a sufficiently large disorder can support 
the generation of entanglement between two distant flux 
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qubits. The main conclusion of this study is that differ- 
ent topics in the fields of metamaterials and localization, 
usually discussed in the classical or many photon level 
can be realized in the few photon limit inside the field of 
circuit QED. 

The paper is structured as follows. In the Sec. II we 
discuss the scattering through a single JJ in the linear 
regime, computing its refiection and transmission coeffi- 
cients. Using these results. Sec. Ill develops the theory 
of transmission lines with periodically embedded Joseph- 
son junctions. We show how to compute and engineer the 
band structure of these photonic crystals and, as applica- 
tion we discuss the implementation of a negative index of 
refraction in two-dimensional arrangements. In Sec. IV 
we study the coupling between qubits and those struc- 
tured lines. We develop an analytical theory that models 
the interaction and dissipation of superconducting qubits 
in a network of J J and transmission lies, within the mas- 
ter equation formalism. This theory is then applied to 
the study of the steady entanglement between two sep- 
arated qubits that interact with a structured line where 
it has been induced disorder. We finish with the conclu- 
sions, while some technical aspects are elaborated in the 
appendices. 



II. JOSEPHSON JUNCTION AS A SCATTERER 

JJs are the most versatile nonlinear element in circuit 
QED. Either alone, or in connection with extra capaci- 
tors or junctions, they form all types of superconduct- 
ing qubits to date^^. Moreover, in recent years they 
have also been used inside cavities to shape and control 
confined photons, dynamically tuning the mode struc- 
ture^^'^^, enhancing the light-matter coupling^^'^^, or 
exploiting their nonlinearity in resonators^^. Junctions 
have also been sugested as control elements for propa- 
gating photons in two different ways. One approach con- 
sist of SQUIDs or charge qubit arrays to control the pho- 
ton dispersion relation forming one dimensional quantum 
metamaterials^ . The other alternative relies on the 
single photon scattering by superconducting qubits^^'^^, 
using the fact that two level systems act as perfect mir- 
ror whenever the incident photon frequency and the qubit 
splitting equals. 

In the following we combine these ideas, providing both 
a uniform theoretical framework to study the interaction 
of Josephson junctions with propagating photons and a 
scalable architecture to construct quantum metamaterials 
by periodic arrangements of these junctions. Just like in 
the case of qubits, we expect that a J J in an open line may 
act as a perfect scatterer of propagating photons, where 
now the resonant frequency is given by the JJ plasma 
frequency. 

In our study we will adopt a bottom-up approach 
starting from the scattering problem of a single junc- 
tion [Fig. 1] that interacts with incoming and outgoing 
microwave packets. The Lagragian for this system com- 



bines the one-dimensional field theory for a transmission 
line with the capacitively-shunted-j unction model for the 
junction^'^'^^'^^ 
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The field ^(x, t) represents flux on the line. The line ca- 
pacitance and inductance per unit length, cq and /q, are 
assumed uniform for simplicity (See Ref. 31 for general- 
izations). The junction, placed at x = 0, is characterized 
by a capacitance Cj and a critical current Ic together 
with the gauge invariant phase, (p 
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where AO is the superconducting phase difference and 
A(r, t) is the vector potential. 

It is convenient to introduce the field (j){x,t) as the 
variations over the static fiux (b^^^x), 
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and a fiux variable S(j){t) associated to the time fiuctua- 
tions for the fiux across the junction S(j){t) := ^(0+,t) — 
^(0_, t) defined as 
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Here (/?^^^ stands for the equilibrium solution for the phase 
and V = (<I>o/27r)(^ = Scj) is the expected voltage-fiux 
relation^^ . 

The fields to the left and to the right of the junction are 
matched using current conservation, which states that 
^left, ^junction and /right are equal at x = 0, 



<I>o 

dx(l){0±,t) = ^Cjif^Icsimp (5) 



These two equations may be formally solved, but the re- 
sult is a complicated nonlinear scattering problem. In 
order to get some analytical understanding of the junc- 
tion as a scatterer, and since we are mostly interested in 
the few photon regime, we will linearize equations (5) 
assuming small fiuctuations in the junction phase d(j)^ 
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with Lj = ^o/{27tIc cos {(p^^"^)). Besides the static fields 
are given by l/lodx(t)''^\x) = /c sin((/p^^^). 
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FIG. 1: (color online) (a) An open transmission line inter- 
rupted by a Josephson junction, (b) Reflection, r, transmis- 
sion, t, and phase of the transmitted beam, Lp = argt, vs. 
incoming photon frequency, in units of the plasma frequency 
ujp. We use Zq/Zj 10. 



In the linearized theory, the stationary scattering so- 
lutions can be written as a combination of incident, re- 
flected and transmitted plane waves: 



i{kx—ujt) _|_ rpQ—'i{kx-\-ujt) 
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(7) 



where is some arbitrary field amplitude, and r and 
t are the reflexion and transmission coefficients, respec- 
tively. We further assume the scattered waves follow a 
linear dispersion relation, cu = vk, which is the same 
outside the junction. Building on the ansatz (7) the co- 
efficients are computed yielding. 
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with the rescaled photon frequency, oo = oo/oOp, with 
ujp = Xj \JLjCj and the impedances of the line and the 
junction, Zj = y^LjJCj and Zq = ^/IqJcq. This for- 
mula, which is analogous to the one for a qubit^^'^^, ex- 
hibits perfect reflection when the photon is on resonance 
with the junction, = cj^, accompanied by the usual 
phase jump across it (cf. Fig. 1). 



III. QUANTUM CIRCUIT CRYSTALS 

We can scale up the previous results, studying peri- 
odic arrangements of junctions both in one and two di- 
mensions. These and other setups^^'^^ can be seen as 
a generalization of photonic crystals to the quantum mi- 
crowave regime, with similar capabilities for controlling 



the propagation of photons: engineered dispersion rela- 
tions, gaps of forbidden frequencies, localized modes, ad- 
justable group velocities'^ and index of refraction, and 
control of the emission and absorption of embedded ar- 
tificial atoms (i.e. improved cavities) 

In the following, we will use the linearized scatter- 
ing theory discussed so far. Regarding possible non- 
linear corrections, we expand the cosine term in (1), 
l/{2Lj)5(^^{l - (27r)V(12/^2)^^2)^ containing both the 
linear contribution and the first nonlinear correction. 
Since we are working in the single-photon regime, we can 
replace 5(j)^ by its fluctuations on the vacuum which are 
proportional to the discontinuity of the wavefunction on 
that point dcj) ^ ^/hZo/2. All together implies a correc- 
tion of 0.2% compared to the linear contribution. If we 
simply view these nonlinear corrections as an inductance 
dispersion we conclude that we can safely neglect them, 
since as we will show in section IV B such a dispersion 
hardly affects the transport properties. 



A. One dimensional circuit crystals 

The simplest possible instance of a quantum circuit 
crystal consists of a unit cell with N junctions that repeat 
periodically in a one dimensional line. The Lagrangian 
is a generalization of Eq. (1), combining the junctions 
together with the intermediate line fields. In the ID case 
there are no additional constraints on the flux and at 
equilibrium cp^^^ = (j)^^^ [x) = minimizes the energy [see 
below Eq. (12)]. The scattering problem is translation- 
ally invariant and its eigensolutions are determined by 
the transfer matrix of the unit cell, Tceii, which relates 
the field at both sides, 4^l,r{x) = cll^rc'^^^ + ^L,i?e~*^^, 
through 
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For a setup with junctions and free lines, the transfer 
matrix has the form Tceii = YliLi^i^i^ where Ti is the 
transfer matrix of the i-th junction and Di is the free 
propagator through a distance di^ 
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The stationary states are given by Bloch waves, which are 
eigenstates of the displacement operator between equiv- 
alent sites. Since this operator is unitary, the eigenvalue 
can only be a phase, ^(xj+i) = exp(zp)^(xj), which we 
associate with the quasimomentum p = kd with d the in- 
ter cell distance. Moreover, as any two equivalent points 
in the lattice are related by the transfer matrix and some 
free propagators, the result is an homogeneous system 
of linear equations whose solution is found by imposing 
det[Tceii(c^)-e^^] =0or. 



2cos(p) = Tr[fceiiM]. 
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FIG. 2: (color online) Photonic crystals with one (a) or two 

(c) junctions per unit cell and their respective energy bands 
(b and d) vs. quasimomentum p. We use Zq/Zj = 10 and 
d — 0.1 A J with A J the typical wave lenght (A = 2'kv/ujj. In 

(d) uj'p/ujp = 0.6 (blue) or 0.9 (gray) and and the distance 
inside the unit cell is O.Ol/d. Notice in (d) that in the lower 
band the gray and blue line are indistinguishable. 



As an example, Fig. 2 shows the dispersion relation 
00 {p) for two simple arrangements. The first one is a line 
with identical Josephson frequency uOp and impedance 
Zj, evenly spaced a distance d [Fig. 2a]. The second one 
is also periodic, but the unit cell contains two junctions 
with different properties, {ujp^Zj) and (a;^,Zj), which 
are spread with two different spacings [Fig. 2c]. We find 
one band gap around u = Up in the first setup, and two 
band gaps around u = Up and uj = ujp in the second, 
more complex case. 

These one dimensional microwave photonic crystals 
have a variety of applications^^. The first one is the sup- 
pression of spontaneous emission from qubits, which is 
achieved by tuning their frequency to lay exactly in the 
middle of a band gap. Another application is the dynam- 
ical control of group velocities. While the width of the 
band gaps is more directly related to the values of Zj 
and the separation among scatterers, their position de- 
pends on the scatterer frequency, ujp. Replacing the JJs 
with SQUIDs^^'^^'^^, it becomes possible to dynamically 
tune the slopes of the energy bands, changing from large 
group velocities (large slope) to almost flat bands (cf. 
Fig. 2d) where photons may be effectively frozen^. Flat 
bands may themselves be used to create quantum mem- 
ories and also to induce a tight-binding model on the 
photons, in the spirit of coupled-cavity systems^^'^^. A 
third application is the engineering of dissipation where 
photonic crystals provide a new arena for theoretical and 



experimental studies. We will focus on this point in the 
last section, studying the relation between disorder, lo- 
calization and entanglement generation in ID quantum 
circuit crystals. 



B. Two-dimensional circuit crystals 

The evolution from one dimensional arrangements to 
two-dimensional or quasi-2D circuit crystals demands a 
careful analysis. The reason for this extra complication 
is that, unlike in ID or tree configurations, phase quanti- 
zation along closed paths introduces new contraints that 
prevent us from gauging away the static phases (/?^^^ and 
fluxes (j)^^\x.,y) in absence on travelling photons. More 
precisely, for any closed path C on the lattice we have 
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Where (fij are the phase differences along each branch 
and <l>ext + ^ind IS the sum of external and induced fluxes 
enclosed by C. The presence of these fluxes may forbid 
an equilibrium condition with all phases equal to zero. 

The physics of our two-dimensional crystals is inti- 
mately related to that of 2D Josephson junction arrays 
(JJAs), a system whose equilibrium and non-equilibrium 
properties have been thoroughly studied in the last 
twenty years^^^^. In particular, we know that JJAs con- 
stitute a physical realizations or the classical frustrated 
XY model, where frustration is similarly induced by the 
fluxes threaded through the 2D plaquettes. 

A proper study of the photonic excitations must there- 
fore begin by studying the static state on top of which 
they will propagate. For that we may rely on the classical 
nonlinear expression for the circuit's energy, built from 
capacitive and inductive terms, where the latter contain 
both the junctions and the (adimensional) mutual induc- 
tance matrix, Lij^ki^^ 



V = -Ej 



ij,kl 



(13) 

The optimization of this problem is a formidable task: 
minimization of (13) subjected at (12) when the induced 
flux are related to the current (phases) through the in- 
ductance matrix. In fact, there is no known solution 
if any DC field is applied. However, let us focus on a 
setup without external fields, then (/p-^^ = 0. This is 
stable against small perturbations and against the quan- 
tum fluctuations induced both by the capacitive terms 
and the travelling photons because the phases are linear 
in the applied field at small fields^^ and they enter on 
second order in the scattering equations [Cf. Eq. (6)]. 



Starting from such stable solution (p[ 



(0) 



0, we can redo the linear scattering theory, which now 
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FIG. 3: (color online) Scattering variables in the two- 
dimensional square lattice. Not shown are the intermediate 
variables: just after the junctions. In Appendix A are denoted 
with a bar. 



contains horizontally and vertically propagating fields 
[Fig. 3], 



(14) 



with k = (kx-,ky). Pretty much like in the one dimensinal 
case, invoking periodicity the solutions are Bloch waves, 




• u=(i,j) (15) 



To obtain the condition for the quasimomentum k = pc/ 
we relate the fields in (i, j) with the ones at {i — 1, j) and 
(z, j — 1) as marked in figure 3. Together with (15) we end 
up in a homogeneus set of linear equations, see Appendix 
A for the explicit calculation . When both the horizontal 
and vertical branches are equivalent, this simplifies to [cf. 
Eq. (11)] 



cos(p^) + cosipy) = Tr[fceii(w)] 



(16) 



based on the transfer matrix along one horizontal or ver- 
tical branch, Tceii(co'), from each elementary plaquette. 



C. Two-dimensional arrays and negative index of 
refraction 



Once we have the possibility of building two- 
dimensional circuit crystals, we can also study the prop- 
agation of microwaves on extended metamaterials, or on 
the interface between them, with effects such as evanes- 
cent waves (i.e. localized modes) and refraction. 

The setup we have in mind is sketched in Fig. 4a, where 
we draw a two-dimensional array of lines with an inter- 
face separating a region with junctions (N) from a region 
where photons propagate freely (F). We may study how 
an incoming wave that travels against the boundary en- 
ters the N region, inducing reflections, changes of direc- 
tion and attenuation. For simplicity we will assume that 



the free region is associated to a vacuum with linear dis- 
persion relation uj'^ = c\{kl. + A:^), where ci the effective 
velocity of light. The region with junctions, on the other 
hand, has an engineered dispersion relation ujN{kx^ky), 
as discussed above. 

When a wave hits the interface between both regions 
it may reflect and refract. The wave that penetrates N 
has to satisfy two constraints: the frequency of the pho- 
tons must be the same as in F and the component of 
the wavevector which is parallel to the interface {py in 
Fig. 4a) also has to be conserved. Both constraints arise 
from a trivial matching of the time {uj) and spatial (py) 
dependence of both waves. Following Ref. 46 both con- 
straints may be solved by inspecting the dispersion re- 
lation in a contour plot [Fig. 4b]. Once the matching 
values of the momenta are found [red in Fig. 4] the ef- 
fective group velocities may be computed to determine 
the trajectory of light. As shown in that plot, in re- 
gions where the dispersion relation is convex, the velocity 
may change orientation and give rise to a refracted angle 
6r = arctan(tanpi^2y cotp2,2y) whose associated index of 
refraction is negative. 

The previous phenomenology has also been proposed 
for a related platform that consists on a two-dimensional 
array of coupled atom-cavity systems. Working in the 
single polar iton subspace, it is possible to derive the dis- 
persion relation for those artificial photons^^ and model 
the array as an effective photonic material. For a sim- 
ilar band structure and interface as the one shown in 
Fig. 4, the effective electrical and magnetic permitivity 
become negative, and one obtains again a negative ref- 
fraction angle^^. The engineering of these counterintu- 
itive refraction processes is of great interest in the field of 
linear optics, as negative indices allow designing perfect 
lenses^^, but the propagation of photons in these mixed 
materials may be interesting also for engineering the dy- 
namics of photon wavepackets, photon routing and 100% 
efficient qubit-qubit interactions — based on the perfect 
refocusing properties of these metamaterials. 



IV. QUBIT-CRYSTAL INTERACTION: 
QUANTUM MASTER EQUATION APPROACH 

So far we have discussed lines with junctions for tai- 
loring photonic transport. In this section we study the 
interaction of these metamaterials with superconducting 
qubits. Modifying light-matter interaction is a corner- 
stone in Quantum Optics. One of the most famous ex- 
amples is the Purcell effect. Confined field enhances or 
dismisses the spontaneous emission for a quantum emit- 
ter. Confinement is usually accomplished by reducing 
the field to one dimensional waveguides or within cavi- 
ties or resonators, see e.g. Ref. 49. Related to this is the 
suppression of spontaneus emission when the transition 
frequency for the qubit is placed inside the gap of a pho- 
tonic crystal^^. While the first is at the heart of current 
circuit QED experiments, the second can be observed 
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FIG. 4: (color online) Two-dimensional circuit crystal, (a) 
Interface between a junction free region (F) and an engi- 
neered bandgap region (N) containing junctions with ujj = 
1.1, Zo/Zj = 0.8 and d = 0.1. Like in the case of polaritons 
in Ref . 46 the square lattice must be tt/ 4-rotated to have neg- 
ative refraction, (b) Dispersion relations for the region N. In 
order to find out the refraction angle for wave that propagates 
from F to N we have to match, on each side of the interface, 
both the photon frequency and the projection of the wave 
vector along the boundary, vi • = V2 • . For large enough 
momenta, (b) shows that the wave gets reversed, that is, the 
index of refraction is negative. 



the Green function for the line, G(x,y,a;), as 

^2 



(19) 



where v = 1/ ^/IqCq is the light velocity and g is the cou- 
pling in a A/2 superconducting resonator with fundamen- 
tal frequency u. This is a very convenient way of express- 
ing the qubit-photon interaction because of two reasons. 
First of all, calculation is simplified to the computation 
of Green functions, wich in our case is particularly easy, 
as it can be derived from the the transfer matrix (9) as 
explained e.g. in Ref. 51 [See Appendix^ for further de- 
tails] . Second and equally important, the strength of the 
coupling is parameterized by a simple number, g^ which 
corresponds to a measurable quantity in qubit-cavity ex- 
periments — ranging from few to hundreds MHz, from 
strong to ultrastrong coupling regimes. 



A. Qubits Quantum Master Equation 

Let us introduce a master equation for Nq identical 
qubits placed in the line at positions xj. The qubits do 
not directly interact, but they will do through the line. 
The qubit Hamiltonian reads, 



(20) 



Interested as we are here in the qubit dynamics, we can 
trace out the transmission line. Assuming, for simplicity 
that the weak coupling qubit-line limit holds, one ends up 
with a master equation for the two qubit reduced density 
matrix^^'^^ 



with our proposal for engineering band gaps. In this sec- 
tion we modelize such light-matter interaction. In the 
weak coupling limit we work a quantum master equation 
and write it in terms of the Green function for the line. 
The latter can be calculated by knowing the scattering 
matrix (10). Finally we give a novel application, as the 
entanglement generation through disorder. To simplify 
the discussion, and without loss of generality, we focus in 
the one dimensional case. 

Let us write the qubit-line hamiltonian^^' , 



line 



Hu 



(17) 



i7q is the qubit Hamiltonian and i^^ne = / duujalja^ the 
line expressed in second quantization. The qubit-photons 
interaction is given by 



iJint = ha"" J dujg{x,uj){a^ + al) 



(18) 



where g{x^uj) is the coupling per mode. In Appendix B 
we show that this coupling can be expressed in terms of 



Oq i 



Here, 



(21) 



(22) 



is the coherent coupling mediated by the line, the so 
called Lamb shift with, 



Ji 



9 



Finally the rates read, 

.(%)' 



7ij = 2tt ImG(ri, r^, e) + A^^ 



(23) 



(24) 



with A the phenomenological non-radiative rate, coming 
from the intrinsic losses of the qubits, see Appendix B 
for details. 
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The simplest situation that is described by this model 
is that of an open transmission line, with no intermediate 
scatterers. In this case the line gives rise to both both a 
coherent and incoherent coupling, quantified by Jij and 
^ij respectively, which depend on the wavelength of the 
photons, the qubit separation and their energies. In this 
case without junctions, G{xi^Xj^uj) = i{v/2uj)e'^^^'^^^^~^^^ 
and thus 

Jij = {v/2uj) cos[ijo/v{xi — Xj)], (25) 
Jij = {v/2uj) sm[ijo/v{xi — Xj)]. 

Note how each of these couplings can be set indepen- 
dently set to zero. This has been used for modify the 
qubits emission going from superradiance and subradi- 
ance^^, and it describes recent results in multi-qubit pho- 
ton scattering^^. 

B. Entanglement through disorder 

With all this theory at hand we move to study a con- 
crete example where we put together structured lines, 
qubits, disorder and entanglement. So far we discussed 
regular (periodic) arr augments of junctions producing an 
ideal photonic crystal. It is feasible to produce them, 
despite fabrication errors, and the junctions within the 
same sample are very similar. Nevertheless it can be in- 
teresting to induce disorder in the scattering elements, 
either statically, intervening in the design or deposition 
processes, or dynamically, replacing the junctions with 
SQUIDs and dynamically tuning their frequencies. Dis- 
order may have a dramatic infiuence in the transport 
properties of the photonic crystal^^. On the one hand, 
the transmission coefficient averaged over an ensemble of 
random scatterers (T) decays exponentially with increas- 
ing length L of the disordered media, similar to Ander- 
son's localization^^. On the other hand disorder fights 
against the interference phenomena that gives rise to the 
existence of band gaps. The consequence of this competi- 
tion will be that a sufficiently large disorder could restore 
the transmission in the frequency range that was origi- 
nally forbidden^^'^^'^^. In what follows we exploit this 
phenomenon in connection to a purely quantum effect: 
the entanglement generation through disordered media. 

Our model setup consists on two well separated ffux 
qubits, Nq = 2 in (21) and (22), which are coupled by a 
quantum circuit disordedered media [Fig. 5] . The qubits 
will be at their degeneracy points and one of them is 
driven by an external resonant classical field: cOd = e, 

= |('^? + '^D + /(e-'"^*<T+ + H.c.) . (26) 

The line, seen now as a quantum bath, has been in- 
terrupted by a set of JJ forming the disordered media. 
The line itself follows our previous scattering theory with 
uniform disorder S in the frequency ujp ^ ujp{l -\- 5) and 
impedance Zj Zj{l-\-5). We use the master equation 



(21) and compute the coefficients in the case of two qubits 
separated by disorder as in figure 5. As demonstrated in 
Appendix C the final expressions read, 

Ji2 = 7oIm(Texp(-zA:2D))/2 , (27) 

accounting for the coherent coupling with the cross- 
dissipation rate 

712 = 7oKe(Texp(-zfc2L))) . (28) 

These interactions compete with the individual decay 
rates of the qubits, 

7ii = 7o(l - Re{ReM-^kD))) + A , (29) 

which includes a phenomenological non-radiative decay 
channel, A coming from the intrinsic losses of the qubits 
[Cf. Eq. (24)]. In all these formulas appear the effective 
rate 70 = TrJi^g^ /vk, the total transmission and reffection 
T and R at the boundaries of the disordered part, and the 
qubit-disorder separation, D. In figure 5 this dependence 
dissapears since the results are drawn at the distance D 
that maximizes the concurrence. 

The physical picture that results is intuitively appeal- 
ing: for the qubits to be entangled, the noisy environment 
should be able to transmit photons, T ^ 0, as both the 
coherent and incoherent couplings depend on it. More- 
over, all photons which are not transmitted but refiected 
add up to the ordinary spontaneous emission rates of the 
qubits, jii. And finally, for a wide parameter range the 
two qubits are entangled also at t ^ 00, in the stationary 
state of the combined system, 9^ ^stationary = 0. We have 
quantified the asymptotic amount of entanglement using 
the concurrence, C, for a variety of disorder intensities in 
a medium which is composed of = 20 junctions which 
are uniformly spread over a distance L = 2A. Fig. 5 
shows the result of averaging 500 realizations of disor- 
der and contains the two ingredients stated above. We 
observe that for zero or little disorder entanglement be- 
comes zero at the band gap, uj/ujp = 1, where photons are 
forbidden due to interference. However, as we increase 
disorder the gap vanishes and entanglement enters the 
region around it. Outside the gap the effect is the op- 
posite: disorder reduces the amount of entanglement, as 
it hinders the transmission of photons. To understand 
the modulations of the plot one must simply realize that 
the value of C mostly depends on the ratio between 712 
and 7ii^^, and these are complex functions of T and i?, 
respectively^^. 

V. CONCLUSIONS AND OUTLOOK 

In this work we have developed an architecture for 
quantum metamaterials based the scattering of travelling 
photons through Josephson junctions. We have shown 
that a single junction acts as a perfect mirror for pho- 
tons that resonate with its plasma frequency. Using the 



8 



/ II I?/sorder~~[J~ 



1 



1 




FIG. 5: (color online) (a) Two qubits connected by a noisy 
environment, (b) Concurrence between the qubits for model 
(21) as a function of frequency and fabrication error (5). We 
simulated a setup with 20 junctions regularly spaced over a 
distance L — 2A, averaging over 500 realizations. We use the 
parameters, Zo/Zj = 10, e = cjd, A = O.470 and / = O.I70. 



scaterring matrix formalism, we have studied the band 
structure of networks of transmission lines with embed- 
ded junctions. We demonstrate that these setups behave 
as quantum metamaterials that can be used to control 
the propagation of individual photons. This opens the 
door to the usual applications of classical metamaterials^ 
such as cloaking or subwavelength precision lenses. In 
particular, as an illustration of the formalism for two- 
dimensional networks, we discussed the observation of a 
negative index of refraction. 

We want to remark that the utility of junction quan- 
tum metamaterials extends beyond the classical regime, 
with interesting applications in the fields of quantum in- 
formation and quantum circuits. Replacing individual 
junctions with tunable SQUIDs opens the door to the 
dynamical control of band gaps, or the generation of flat 
bands, which is useful for stoping light, implementing 
quantum memories and what would be the equivalent of 
coupled cavities arrays. 

Two important applications of this tunability are engi- 
neering of disorder and dissipation. In the first case the 
focus is on the photons that travel through the network, 
while in the second case the focus is on how this network 
acts on few-level systems that are embedded in them. 
We combine both approaches by developing the theory 
for multi-qubit interactions in a quantum metamaterial. 
The resulting master-equation formalism combines the 
effects of spontaneous emission in the artificial material, 
with the interaction mediated by the exchange of pho- 
tons. We show that two competing effects — Anderson 
localization suppresses transport, but disorder populates 
the band gaps with localized states — , lead to the gener- 



ation of stationary entanglement in these setups. 

We strongly believe that this architecture is within 
reach for the experimental state of the art. Building on 
very simple components, it offers a great potential both 
for quantum information with flying microwave qubits 
(photons), and for the static and dynamic control of sta- 
tionary qubits. In the near future we wish to explore 
the application of this technology as a replacement for 
the coupled-cavity architecture, where the one- or two- 
dimensional network replaces the cavities, offering new 
possibilities of tunability and variable geometry. 

Finally, we want to remark that shortly after the 
rewrite of this manuscript, a related work appeared that 
develops a similar formalism for Josephson junctions em- 
bedded in transmission lines^^. 
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Appendix A: Transfer Matrix in the 2D case 

We detail here the calculations needed to obtain con- 
dition (16) in the main text. The idea is to relate the 
horizontal and vertical fields (14) on both sides of each 
junction, as shown in Fig. 3. Introducing the vectors 



(h,v) _ / 



(Al) 



the fields at both sides of the junction are related, see 
Eq. (9), 

^^=T^(co)^f!},^, (A2) 
^§^=T^(oj)^^l„ (A3) 

Finally we get the final w-^'^^ by resorting to continuity 
and current conservation in the corner. It is convenient 
then to define the vectors e!^ = (1, ±1). Using this nota- 
tion, the continuity condition is written as follows 



t rpih)(h) _ t rpiv)(v) _ t (h) _ t (V) 



T^"^M<r^"^ - T 



and current conservation at the corners reads 



(A4) 



etT^^,^ + ^t!^^1, = e*_w(f + e*_wg). (A5) 
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writting now w-^'^^ as Bloch waves (15) we end up 4- 
coupled homogeneus set of linear equations: 



with, 



l{u;) 



t ih) 



(A6) 



/ 1 

1 



ika:„t rp(h) 





3* 



-ik 



y 1 _ e— xet_Tv-e+ 1 
together with the relations due to the scattering matrix (10) properties, 



— e 
1 - e 







-iky^t rpiv) 
^-^ceir 



r 



(A7) 



(A8) 
(A9) 



Putting altogether we have that det[M] = yields 
the generalized condition for the two-dimensional case. 
In the simplest case of fully symmetric configuration: 



"-cell 



:= Tceii we simply have [Cf. Eq. (11)], 
cos(p^) + cos(p^) = Tr[Tceii(cj)] (AlO) 



Appendix B: Modelling qubit-line interaction and 
Master Equation 

In this appendix we develop the model for the qubit- 
line interaction. We will focus on fiux qubits for the 
sake of concreteness, but the results are analogous for 
other qubits. Besides we will discuss the master equation 
governing the qubits dynamics. Finally we rewrite the 
formulas in terms of the Green function. 

For fiux qubits the coupling is inductive and can be 
written in circuit and/or magnetic language as. 



ML 



qubit 



X /, 



line 



(Bl) 



here M stands for the mutual inductance, /qubit and /une 
are the currents and fi is the magnetic qubit dipole, while 
B is the magnetic field generated in the cavity. 

The current in the line is given by /une = i^dx(l){x). We 
will expand this field using normal modes, Un{x)^ follow- 
ing the usual quantization (j){x^t) = ^Uk{x)qk{t)^ but 
imposing that Uk are dimensionless^^ and satisfy the or- 
thonormality condition / coUk uidx = CrSki with the av- 
erage capacitance Cr := / codx. Expressing the canoni- 
cal operator qk in the Fock basis Qk = {a\^ak)\/h/2ujriCr 
gives us the final expression 



^line 



=-y 



dxUk{x){al + Uk) . 



(B2) 



The magnetic field-current relation is given by ^une = 
Mo/iine/7r(i, with d the distance between plates in the 
coplanar wave-guide. The quantized magnetic dipole for 
the qubit can also be expressed in terms of the qubit area, 
A and the stationary current Ip as /i = IpAa^ . Putting 
alltogether we find the interaction Hamiltonian (Bl) 

Hint = IpA^^^ da, yZ^=^^^k{x){al^ak) . 

7ra to — y^k 



(B3) 

We can introduce the coupling strenght per mode with 
frequency cjq^^, 



hg IpA 



Mo 



1 



(B4) 



Grouping the constants and using the expression for 
Eq. (B4), we rewrite (B3), 



Hint = h^v^' 



2L 



^=dxUk{x){al + a/e) 



(B5) 

where v = I/V/qCq the light velocity in the line and ujq 
the fundamental frequency of a cavity with a given g and 
L is the Lenght. As expected the above expression is 
nothing but the spin boson model. 



1. Green Function formalism 

It turns out useful to rewrite (B5) in terms of the Green 
Function for the line. We begin the discussion by re- 
calling the field wave equation. Equivalently to layered 
photonic crystals, it is sufficient to work the case of ho- 
mogeneous line, since the problem we are dealing with is 
piecewise homogeneous^^. For fiux qubits the coupling is 
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through the hne current [cf. Eq. (Bl)]. Thus it is more 
convenient to discuss the wave equation for the mode 
derivatives 



(B6) 



with the orthogonality condition, 



/ 



dxdxUkdxUk' = Lk'^Skk'- (B7) 



The Green function for this Sturm-Liouville problem 
reads 

dlG{x, x\ u) + -;^G{x, x\ u) = -S{x - x'). (B8) 

It is pivotal the relation [cf. Eq. (8.114) in Ref. 63] 

IT ^ dxUk{x)dxul{x') 



ImG(x,x\cj) = 



L 2 



5{u - Wfe) . 

(B9) 



By rewriting (B5) in the continuum limit, 

i^int = her'' j (lujg{x,uj){a^ + a^) (BIO) 
with g{x^uj) [combining (B9), (B5) and (BIO)] 



\g{x,uj)\^ = 27T^lmG{x,x,uj) , (Bll) 
i.e. equations (18), (19) in the main text. 



2. Quantum Master Equation 

Setting the temperature to zero (typical experiments 
are at the mK while frequencies are GHz) the qubit dy- 
namics, after integrating the bosonic modes, is given by 
the standard master equation in Linblad form wich as- 
sumes weak coupling between the line and the qubit^^'^^. 



dtp 



: [^qubit + ^LS , p] 



(B12) 



where {, } is the anticonmutator. In the equation we have 
distinguished the contribution to the decay rates com- 
ing from the qubit-line coupling, F^j, from other noise 
sources affecting the qubits, denoted with a phenomeno- 
logical strenght A. The explicit expressions for the Ti^j 
are e.^.^^'^^ 



7(e)l' 



(B13) 



X = x/ X = 



X = L 



FIG. 6: Sketch for the Green function calculation. The "black 
box" is characterized by transmission and reflection coeffi- 
cientes. The source (Dirac Delta) is represented by the ring 
attached to the line. 



Finally, we also have to consider the Lamb Shift 



with 



V 



vHm\giu)\'^ 



(B14) 



(B15) 



where V[ ] means principal value integral. 
By defining. 



(B16) 



together with (Bll) we get coefficients (23) and (24) in 
the main text. 



Appendix C: Green Function for an arrangement of 
scatterers 

In the following we find G{xi^Xj^uj) for the problem 
discussed in the main text. We show that G{xi^Xj^uj) 
is written in terms of reflection and transmission coef- 
ficients, R and T respectively. We use this to express 
the decays and cross couplings, jij and J12, in terms of 
the scattering parameters, making explicit the connec- 
tion between the photonic transport in the line and the 
dynamics for the qubits coupled to it. 

In our case, two qubits placed at possitions Xi and 
X2 with a set of junctions in between (see figure 6) 
G{xi^Xj^uj) can be computed as follows. The equation 
for the Green function (B8) is a field equation with a 
source (because of the Dirac delta) at x = x^ The junc- 
tions cover the region from x = to x = L, therefore 
xi < and X2 > L. This situation is analogous to have 
a boundary with reflection R and transmission T, as de- 
spicted in figure 6. In this situation the Green Function 
is given by [Eqs. (2.34) and (2.35) in^^]. 



G(x, x\ uj) 



2k 



2k 



-ik{x-x') _ j^Q-ik{x-\-x')^ X < x' 
^ik{x-x') _ Ji^-ik{x+x')^^ x' < X < 

x> L 



i rp ik(x — (x' -\-L)) 



We remind that the minus sign in front of the R above 
comes because the Green function in (B9) is given in 
terms of dxUk. 
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Finally, the coefficients in the master equation read 
[Cf. Eqs. (24) and (23)] 

rjj(a;qubit) = 27r^^^ImG(r^-,r^-,a;qbuit) 



ri2(wqubit) = 27r^^^ImG'(ri,r2,Wqbuit) 

= 27r^^^Re(re-*'=(^i-^^)) (C2) 

V Zrv 

The last obstacle to write the master equation is to 
perform the integral in (23). Here we made use of the so 



called Generalized Kramers-Kroning relation^"^, namely 



V 



f 

Jo 



uj^ ImG(rj, r/e, cj) 



Thus, 



Ji2 = ^^^^Im(Te-^^("^-"^)) 



V 2k 

Introducing the definition 

70 TT- 



.(%)!i 

V k 



|-^jf^ReG(rj,rfe,w) 
(C3) 

(C4) 
(C5) 



we end up with the expressions used in the main text. 
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